Orbital Polarization in Itinerant Magnets 
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We propose a parameter-free scheme of calculation of the orbital polarization (OP) in metals, 
which starts with the strong-coupling limit for the screened Coulomb interactions in the random- 
phase approximation (RPA). For itinerant magnets, RPA can be further improved by restoring the 
spin polarization of the local-spin-density approximation (LSDA) through the local-field corrections. 
The OP is then computed in the static GW approach, which systematically improves the orbital 
magnetization and the magnetic anisotropy energies in transition-metal and actinide compounds. 

PACS numbers: 71.15.Mb, 75.10.Lp, 71.15.Rf, 75.50.-y 



An electron in solid can carry spin (Ms) and orbital 
(A/i) magnetic moment. For weakly correlated systems, 
the problem of spin magnetism alone can be formulated 
in the fully itinerant fashion, meaning that the effect 
of other electrons onto a given one can be described 
by an exchange-correlation field (or spin polarization). 
The field is typically evaluated in the model of homoge- 
neous electron gas, in the basis of plane waves, which is 
a limiting case of the extended Bloch waves. This con- 
stitutes the ground of the Kohn-Sham (KS) formalism 
within LSDA 1], which works exceptionally well for the 
magnetic spin properties of many transition-metal and 
actinide compounds. They form an extended group of 
what is currently called the "itinerant electron magnets" . 

The orbital magnetism is an atomic phenomenon. In 
the majority of cases, it is driven by the spin-orbit inter- 
action (SOI), being proportional to the gradient of the 
one-electron potential, VV, which is large only in a small 
core region close to the atomic nucleus. Furthermore, 
the angular momentum operator, L^, does not commute 
with V. Generally, is not an observable quantity, 
except the same core region, where V is spherically sym- 
metrical. Therefore, it is more natural to formulate the 
problem in the basis of Wannier orbitals {(j>a} (a being a 
joint spin-orbital index), localized around each atomic 
sites 2]. Then, the orbital moment ML=TrLs{L^n} 
is specified by the local density matrix n=||riQ^|| (Tr^s 
being the trace over spin and orbital variables), where 
nai3=Yjinidaid}p^, dai={4>a\tpi) is the projection of KS 
eigenstate tpi onto is the KS occupation number 

corresponding to the eigenvalue e^, and the joint index i 
stands for the spin, band, as well as the position of the 
k-point in the first Brillouin zone (BZ). 

In an analogy with the spin polarization for itiner- 
ant magnets, one can think of an OP: an exchange- 
correlation field in KS equations, which couples with M/,. 
Despite a genuine interest to the problem and wide per- 
spectives of their potential applications, the theories of 
OP in metals are still in a developing "semi-empirical" 
stage, as they largely depend on the input parameters, 
which are typically chosen to fit the experimental data. 



Although majority of researches agree that OP is con- 
trolled by intra-atomic interactions, which are strongly 
screened in metals, the details of this screening as well 
as the form of the OP itself remains to be a largely un- 
resolved and disputed problem dH H 0,0,13 ■ 

Therefore, there are two important questions, which 
we would like to address in this work, (i) How the bare 
on-site interaction Uaf3'yS={4>a4''y\^ /''^I2\4'i3(t's) between d- 
or /-electrons is screened in metals? What is the main 
mechanism of this screening? (ii) Is there any simple 
and reliable way to evaluate this screening in ab initio 
calculations of OP? 

In the atomic limit, the full matrix u^Hua/j-yil is con- 
trolled by a small number of Slater integrals {F'^}. Then, 
there is an old empirical rule , which states that in met- 
als, the screening affects mainly which contribute to 
the Coulomb matrix elements Uaa-yy- Other Slater inte- 
grals, which control the exchange and nonsphericity of 
Coulomb interactions do not change so much. 

First, we argue that the same type of screening can 
be naturally obtained in RPA, in the fully deterministic 
fashion. The screened interaction in RPA 0, ^| , 



U{uj) = 1 - uP{uj) 



(1) 



depends on the polarization -P=|| Pa/3715 1|, which is treated 
in the approximation of noninteracting KS quasiparticles: 
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The w-dependence of P contributes mainly to the 
redistribution of the spectral density, whereas the uj- 
integrated ground-state properties are controlled by U = 
C/(0). Therefore, we consider only the static limit, in 
which RPA describes the screening of u caused by the 
relaxation of {tpi} upon removal or addition of an elec- 
tron in terms of the perturbation-theory expansion [l3| . 

The simplest toy model, which illustrates the physics, 
may consist of two spin-polarized bands, formed by (I) 
and zx (2) orbitals. The model is compatible with the 
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orbital magnetization in the (001) direction. Adopting 
the following order of orbitals (within one spin channel) : 
af3 {■^5)— 11, 22, 12, and 21, it is easy to show that 
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(3) 



and 



where u=F^+±F^+§-,F\ j = I^F^+^^f\ 
u'—u~2j. Due to the orthogonality of the yz and zx or- 
bitals, the Coulomb {aP= 11, 22) and exchange {a/3= 12 
21) matrix elements are fully decoupled from each other 
In order to illustrate the main idea of RPA-screenin; 
P can be taken in the form Paii-yS^P^aS^p-y 
which yields: U=[u-{v?-u''^)P]/[{l-uPf-{u 
U'^u' /[{l-uPf-{u'P)\ and J=i/[l-jP]. There is 
certain hierarchy of bare interactions, and for many 
metals the screening of u and j falls in the strong- 
and weak-coupling regime, respectively, so that u|P|3>l 
while j|P|<l 111. This yields: J7~-1/(2P)+2J, 
?7'~— 1/(2P), and J~j. Thus, this is the inverse 
polarization, which plays a role of effective Coulomb 
interaction in metals [l^. U is strongly screened and 
does not depend on the value of hare interaction. On 
the other hand, J is insensitive to the screening. The 
multiplier 1/2 in the expressions for U and U' stands 
for the orbital degeneracy. The result can be easily 
generalized for an arbitrary number of orbitals M {M— 
5 and 7 for d- and /-electrons, respectively), which 
yields C/'~— 1/(MP) In this case, in order to justify 
the strong-coupling regime, it is sufficient to have a 
milder condition, uM|P|^l, which naturally explains 
the empirical rule 0. 

All these trends are clearly seen in the first-principles 
calculations for realistic materials shown in Fig. ^ [l5j , 
where all Slater integrals except P° were calculated in- 
side atomic spheres, and P*^ was treated as a parame- 
ter. When P° increases, the effective interactions quickly 
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FIG. 1: Effective Coulomb and exchange interactions in RPA 
versus bare Slater integral P" for 3d-states in bcc Fe and 5/- 
states in uranium sulfide. The symbols denote the matrix ele- 
ments corresponding to different representations of the point 
group Oh- The calculations have been performed in the fer- 
romagnetic state without spin-orbit coupling. 

reach the asymptotic limit F^—xx, where U is fully de- 



termined by details of the electronic structure, through 
the polarization P, and do not depend on F'^. This re- 
moves the main ambiguity with the choice of interaction 
parameters for metallic compounds. Since P depends on 
the local environment in solid, the screened interactions 
can be different for different types of Wannier orbitals 
(e.g., Cg and t2g for d-electrons in the cubic environment). 

Thus, OP in the itinerant magnets can be naturally 
evaluated in the framework of an universal parameter- 
free scheme based on the strong coupling limit for the 
matrix of effective Coulomb interactions U. The self- 
energy, incorporating the effects of OP, can be calculated 
within static approximation in the GW method |1(T| | : 



7(5 



(4) 



The proper correction to the KS Hamiltonian in LSDA 
is controlled by Afi—h—^j^ X]r=o TrLs{<5'r"-}<5'r- It is ob- 
tained after subtracting the charge (r=0) and spin {r= 1, 
2, and 3) density elements of n, which are already taken 
into account in LSDA (ctq being the unity matrix, and 
(7i , (T2 , and (73 being the Pauli matrices of the dimension 
2M). Therefore, in the actual calculations we uses the 
change of the self-energy AE, which was obtained after 
replacing n by An in Eq. The problem was solved 
self-consistently with respect to An. 

The validity of the strong-coupling approach is well 
justified. So, the effective Coulomb interaction between 
t2g electrons in bcc Fe can be estimated in RPA as 1.50, 
1.47, and 1.37 eV for po= oo, 21 eV (the bare Slater 
integral inside atomic sphere), and 4.5 eV (the value ob- 
tained in the constraint-LSDA, which includes the screen- 
ing by the sp-electrons 0), respectively. Thus, even if 
one takes the lowest estimate P"=4.5 eV, the additional 
approximation F'^—^oo within RPA would overestimate 
U by less than 10%. For /-electrons, this error is even 
smaller due to the higher orbital degeneracy. 

However, this is not the main source of the error. A 
more fundamental problem is related with the RPA it- 
self, which typically underestimates the spin polariza- 
tion A^^^=Tris{StT3}, meaning that even for the up- 
per limit in RPA, corresponding to P'^^oo, the effective 
Coulomb interaction is overscreened and underestmated. 
For example, had we replaced the spin part of LSDA 
by the one of RPA, the spin moment would be underesti- 
mated. Obviously, this would destroy the most attractive 
point of LSDA for itinerant electron magnets. Therefore, 
there is certain inconsistency in the RPA approach. 

RPA can be improved by introducing the local- 
field factor g, which incorporates the effects of 
exchange-correlation hole for the polarization matrix: 
(P)-i=(pR-PAv-i_^. Other corrections can be formally 
reduced to g [ig] • Our goal is to find such a correction to 
the matrix of effective Coulomb interactions, which after 
substitution in Eq. I^J would yield the same spin polar- 
ization as LSDA (A'"^^^). In order to do so, we search 
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g in the form of local diagonal matrix: ga(3^s=g5aj35^s- 
Then, the asymptotic part of the effective Coulomb inter- 
action and the self-energy can be easily recalculated using 
Eqs. (Q) and Q , respectively, and the unknown parame- 
ter g is obtained from the condition Tri5{S(T3}=A'^^^^, 
which is solved self-consistently together with the KS 
equations. In the following, this method will be referred 
to as corrected- RPA (c-RPA). 

Let us consider first the canonical example of ferromag- 
netic transition metals (Fig. EJ, where Ml is small and 
typically regarded as a small perturbation to the spin- 
dependent properties. Ms and Ml can be measured sep- 




FIG. 2: (Color online) Spin (light blue area), orbital (dark 
red area), and total (full hatched area) magnetic moments in 
ferromagnetic transition metals. The experimental data are 
taken from Ref. P|. The inset shows the Fermi surface of fee 
Ni in the c-RPA approach. 

arately using the x-ray magnetic circular dichroism com- 
bined with the spin and orbital sum rules Despite 
an apparent simplicity, LSDA encounters a wide spec- 
trum of problems for bcc Fe, hep Co, and fee Ni. We will 
argue that many of them can be systematically corrected 
by applying consequently RPA and c-RPA techniques. 
LSDA has certain tendency to overestimate Ms in bcc Fe 
and underestimate M^, while RPA and especially c-RPA 
substantially improve the LSDA description and yield a 
good agreement with the experimental data. The values 
of Ms [Ml) obtained in LSDA, RPA, and c-RPA are 2.26 
(0.04), 2.21 (0.05), and 2.20 (0.06) hb, respectively, to be 
compared with the experimental moments of 2.13 (0.08) 
Ms lal- Hep Co has the largest orbital moment among 
pure transition metals (Mi=0.14^B), which is strongly 
underestimated in LSDA (Ml=0.08/l<b)- The situations 
is substantially improved in RPA {Ml=Q.1Qixb) and c- 
RPA {Ml—Q.11^b)- Fee Ni is a rare example of fer- 
romagnetic systems for which 71/^=0. 05/LtB is well re- 
produced already in LSDA. Both RPA and c-RPA pre- 
serve this good feature of LSDA and do not substantially 
change Ml ■ However they do change the electronic struc- 
ture of fee Ni. Namely, the form of Fermi surface (FS) 
of fee Ni has been intensively discussed in the context 
of the magnetocrystalline anisotropy energy (MAE). It 
was argued that the reason why LSDA fails to reproduce 
the correct (111) direction of the magnetization is related 
with the second pocket of the FS around the X-point of 
BZ, which is not seen in the experiment The experi- 
mental FS can be reproduced in the LSDA-)-?7 approach. 



by treating U as an adjustable parameter 0. Therefore, 
it is important that the same problem can be successfully 
resolved both in RPA and c-RPA, without any adjustable 
parameters. The calculated FS, which reveals only one 
pocket around the X-point, is shown in the inset of Fig.|21 
The uranium pnictides (UAT, where X= N, P, As, Sb, 
and Bi) and chalcogenides {X= S, Se, and Te) are ones of 
the most studied actinide compounds. They crystallize in 
the rock-salt structure. All chalcogenides are ferromag- 
nets, whereas the pnictides have type-I antiferromagnetic 
structure, which may also transform into a multi-k struc- 
ture. The basic difference from the transition metals is 
that Ml in actinides, which can be extracted from the 
analysis of magnetic form factors 0, , is very large and 
typically dominates over Ms- According to the third 
Hund rule, Ms and Ml in UX are coupled antiferro- 
magnetically. As the U-U distance increases, the U(5/) 
states become more localized, and all magnetic moments 
increase monotonously from UN to UBi and from US 
to UTe (Fig. OJ. UN and US are usually classified as 
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FIG. 3: (Color online) Magnetic moments in uranium 
pnictides (top) and chalcogenides (bottom). The pnictides 
(chalcogenides) have been computed in the type-I antiferro- 
magnetic (ferromagnetic) structure with (001) ((111)) direc- 
tion of the magnetization. The symbol 'exp' shows the results 
of neutron diffraction, which were separated into spin and 
orbital contributions for US (Ref. (i|) and UAs (Ref. 
Other notations are the same as in Fig.|5| 

itinerant magnets. However, the role of intra-atomic cor- 
relations is expected to increase for the end-series com- 
pounds. Obviously, the real ah initio scheme does not 
know whether the system is itinerant or not. Therefore, 
it is important to test both RPA and c-RPA methods 
for all considered compounds in order to see how they 
will work for the materials with the different character 
of the 5/-electrons. The orbital moments are systemat- 
ically underestimated in LSDA. The error is really large 
so that the total magnetic moments are typically off the 
experimental values by 20-50%. RPA systematically im- 
proves the LSDA description. However, it is not enough, 
and for many uranium compounds it is essential to go 
beyond RPA. For these purposes, c-RPA works excep- 
tionally well and further improves the RPA description. 
Particularly, we note an excellent agreement with the ex- 
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perimental data for X= S, P, and As. For the end-series 
compounds {X= Te, Sb, and Bi) the agreement is not 
so good, signaUing at the necessity of more radical im- 
provements, involving both orbital and spin polarization 
of LSDA. However, even for these complicated systems, 
c-RPA is a big step forward over conventional LSDA. 

Finally, let us discuss applications for the MAE. We 
consider two characteristic examples: CoPt and US. The 
ordered tetragonal CoPt alloys is a promising candi- 
date for magnetic recording applications. An intrigu- 
ing point is that although LSDA underestimates Ml, 
MAE is reproduced surprisingly well (Fig.Q j^. There- 
fore, the "correct" OP in CoPt should affect only M^. 
This requirement is well satisfied both for RPA and c- 
RPA. The orbital moments systematically increase in the 
direction LSDA-^RPA-+c-RPA to reach Mf°=0.14^B 
and A/£'=0.07/iB. The anisotropy of Ml also increases 
(mainly at Co-sites). However, the MAE does not change 
so much because of large cancellation of on-site interac- 
tion energies associated with Co and Pt sites jSl]- 




Pt-;itom 
LSDA KPA ppa 




FIG. 4: (Color online) Magnetocrystalline anisotropy energy 
(Ai5) and the anisotropy of orbital magnetization (AMl). 
For each quantity, the anisotropy is defined as the difference 
between values corresponding to the (100) and (001) (CoPt), 
and (100) and (111) (US) directions of the magnetization. 
The experimental values are taken from Ref. 8] (CoPt, at 
293 K) and Ref. [l^ (US). For CoPt, the values of Ms and 
Ml in the (001) direction are shown in the left part of the 
figure. Other notations are the same as in Fig. |5| 

US has the largest MAE among cubic compounds ^3 ' 
which is underestimated in LSDA. The situation is cor- 
rected in c-RPA, at least qualitatively. R is curious 
that MAE "anticorrelate" with the anisotropy of or- 
bital magnetization, which decreases in the direction 
LSDA— >RPA— >c-RPA. However, this is not surprising, 
because in cubic compounds, MAE is the forth order ef- 
fect with respect to SOL Therefore, there is no simple 
relation between A_E and AM^ and the main correction 
to MAE in c-RPA comes from the change of the on-site 
interaction energy. 

In summary, we have argued that the problem of OP 
in metals can be naturally formulated "from the first 



principles" , by considering the strong-coupling limit for 
the screened Coulomb interactions. In the present work, 
the screened C/ was computed only once: in LSDA and 
without SOI. An important extension would be a self- 
consistent determination of U , which would incorporate 
the effects of OP into the screening, (i) It could improve 
the description of some itinerant actinide compounds 
(e.g., UN) for which the spin polarization is small, and 
the screening is strongly influenced by SOI. (ii) Since the 
OP affects the KS eigenvalues {ci}, which stand in the 
denominator of the polarization matrix ^ , the screening 
is expected to decrease. This could extend the applica- 
bility of the proposed method for materials with more 
localized 5/- and 4/-electrons. 
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